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Spectral analyses are given of the groundwave portion of 33 sferic waveforms recorded 
from cloud-to-ground lightning discharges which occurred at distances ranging between 
about 150 and 600 kilometers from Boulder, Colo. Frequencies of peak energy lie between 
5 and 20 kilocycles per second, which agree favorably with other published results. The 
average value of energy calculated from the groundwave pulses was found to be 26, 600 
joules, which is lower than values derived from other experiments. Various parameters, 
such as the peak amplitude and duration of the first half-cycle, are related to the radiated 
energy of the stroke. 



1. Introduction 

Appleton et al. (1926) were probably the first 
workers to utilize the cathode-ray oscilloscope in 
studying the electromagnetic energy emitted from 
lightning- discharges. Following the development 
of radio communication techniques, it became 
evident that the lightning discharge was the source 
of interference to radio communication circuits. 

In recent years, the radio signals emitted from 
lightning discharges, culled atmospherics, or sferics, 
have been utilized as a source of signals in propaga- 
tion studies at very-low radiofrequencies. 

Various workers in England have been particularly 
active in the utilization of sferics in propagation 
studies, and it would be difficult and exhaustive 
to properly credit all of their accomplishments. 
However, as the work continued, it was learned 
that some atmospheric waveforms could be inter- 
preted in terms of the reception of a series of pulses 
due to successive reflections between the earth 
and the ionosphere. It became evident that the 
character of the atmospheric waveform observed 
at great distances was materially altered by the 
propagation effects, while at shorter ranges, the 
waveform was more representative of the source 
function. 

At the present time the use of atmospherics in 
propagation studies at the very-low radiofrequencies 
is common. NBS Boulder Laboratories have es- 
tablished a network of stations [l] 1 for this purpose, 
and the material presented in this paper resulted 
from the simultaneous observation of atmospherics 
at the Bureau recording sites. 

The objective of this paper is to present data 
on Fourier spectra of the electric field strength of 
signals radiated from return-stroke lightning dis- 
charges. The radiation spectra of lightning dis- 
charges were determined from observations of the 
groundwave portion of the atmospheric. The pre- 
cautions taken in identifying and utilizing ground- 
wave pulses for this purpose are described in this 
paper. 



i Figures in brackets indicate the literature references a1 the end of I his paper. 



2. Equipment and Collection of Data 

Stations at Boulder, Colo.; Salt Lake City, Utah; 
and Palo Alto, Calif.; were equipped to record the 
directions of arrival of atmospheric pulses, the 
pulse waveforms, and timing marks from which the 
times of arrival of individual pulses could be deter- 
mined. Atmospheric waveforms used in the spectral 
analyses were recorded from vertical antennas. The 
overall amplitude response of the vertical antenna 
channel was constant within ±1 db over the hand- 
pass and sloped to 3-d b cutoff points at frequencies 
of 1 and 100 kc. The phase response of this channel 
closely approximated a linear function of frequency 
within the bandpass. 

A pair of vertical electrostatically shielded loop 
antennas, arranged at right angles to each other, 
were used in a direction-finding system at each 
station. The overall bandpass of each loop-antenna 
channel extended from 1 to 100 kc, and closely ap- 
proximated that of the vertical-antenna channel. 
The direction-finding indications were obtained 
using the entire bandwidth of the loop-antenna 
channels. It was determined in earlier experi- 
ments [2] with wideband direction finders, that in 
many cases the groundwave component of the 
atmospheric could be discerned from the skywave 
components by virtue of differences in polarization 
and times of arrival of the groundwave and skywave 
pulses. 2 A time resolution of approximately 10 
Msec is afforded by the wideband direction-finder, 
which is adequate to resolve the groundwave and 
skywave pulses. The direction of arrival of the 
atmospheric, as indicated by the vertically polarized 
groundwave component, should be virtually free 
from error due to the presence of horizontally polar- 
ized components in the atmospheric. 

Timing marks derived from a secondary frequency 
standard, synchronized with the standard timing 
emissions from station WWV, were recorded on the 
atmospheric waveform records. Using these timing 

2 Excellent discussions regarding the location of sferic sources through the use 
of radio direction-finding have appeared in the literature; Horner [3], in particu- 
lar, has discussed the probable sources of error in direction-finding, including 
that resulting from the reception of horizontally polarized components of iono- 
spheric waves. 
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marks, the time of arrival of atmospherics could be 
determined at each station to an accuracy of about 
±1 msec. Atmospheric waveforms resulting from 
a particular lightning discharge were located on the 
photographic records made at each station by virtue 
of their times of reception and directions of arrival. 

The sweep of the oscilloscopes used in recording the 
atmospherics was activated by signals having 
amplitudes of about 50 mv/m or greater. Signals of 
smaller amplitude did not activate the sweep and 
therefore were not recorded. This feature provided 
some discrimination against atmospherics of rela- 
tively low amplitude; such as those of distant origin 
as well as some local cloud-to-cloud and precursory 
discharges. To minimize loss of early detail in the 
waveform, a 24-,usec delay line was used to delay 
recording the signal until after the sweep was initi- 
ated. The dynamic range of the recording equipment, 
in terms of the strength of a vertically polarized wave 
incident at the monopole antenna, extended from 0.2 
to 6 v/m. Most of the analyzed atmospherics had 
peak field amplitudes approaching 1 v/m. 

The locations of lightning discharges were deter- 
mined by triangulation, using the direction of arrival 
indicated at each station. Due to the geometric 
arrangement of the 3-station network, the locations 
of lightning discharges which occurred near a record- 
ing station are relatively inaccurate; however, for 
distances greater than about 150 km, the error is 
estimated to be equal or less than 10 percent of the 
range. 

It was desired to identify waveforms on the photo- 
graphic records that resulted from vertical lightning 
discharges. Two tests that were used in an attempt 
to achieve this selection are described below. The 
relative amplitude of the vertical component of the 
electric field which would be registered at the appro- 
priate ranges from a vertical discharge was calculated 
using range information available from the direction- 
finding "fixes," and from groundwave propagation 
curves in the literature [15]. The source was con- 
sidered to be vertically polarized if the relative 
values of the observed and calculated pulse ampli- 
tudes agreed within about 20 percent. 

If a vertically polarized wave arrives at grazing 
incidence, the direction-finder response will approxi- 
mate a straight line inclined at an angle proportional 
to the azimuthal direction of arrival of the wave. 
An incident wave consisting of both horizontal and 
vertical components of the electric field will cause 
an elliptical direction-finding pattern to be produced. 
The direction-finder responses to the atmospheric 
pulses which were selected for analyses indicated 
that the groundwave pulses were vertically polarized, 
and the first-hop skywave pulses were essentially, 
but not entirely, vertically polarized. Atmospheric 
waveforms which satisfied both the "amplitude" 
and "ellipticity" criteria outlined above were as- 
sumed to result from vertical lightning discharges. 

No sferics were analyzed which originated at 
distances less than 150 km from Boulder, due, in 
part, to limitations in the accuracy of "fixes" at 
short ranges, and also to minimize contributions by 



the induction field at frequencies near 1 kc. It was 
desirable to limit the maximum observation range 
in order to minimize the distortion of the ground- 
wave pulse by propagation and to reduce the inter- 
ference caused in the groundwave by the first-hop 
skywave pulse. 

The waveforms of two sferics, representative of 
those used in this analysis, are reproduced in figure 
la. The range to the source is approximately 165 
km from both sferics 1 and 2. The groundwave 
portion of sferic 1 is believed to comprise the first 
110 //sec of the approximate 450-/xsec time base. 
The earliest evidence of the arrival of the first-hop 
skywave pulse in sferic 1 occurs at 147 yusec. Note 
that the initial half cycle of the skywave pulse is 
negative, or opposite to the polarity of the first 
half -cycle of the groundwave pulse. This reversal 
in polarity of the initial part of the first-hop skywave 
pulses evident in figures la and lb is a consequence 
of the pseud o Brewster angle in ionospheric reflec- 
tion. This phenomenon was discussed in an earlier 
paper [1]. 

The height of the ionospheric reflection was 
calculated to be about 65 km, corresponding to a 
range of 165 km and a 147-jusec time interval between 
the reception of the groundwave and the first-hop sky- 
wave pulses. In this calculation, the ionosphere 
was assumed to be a homogeneous ionized medium 
with a sharp lower boundary [4]. 

Note that sferic 2 is similar in character to sferic 
1, but the groundwave is of larger amplitude and 
longer duration. The first-hop skywave pulse was 
received before the entire groundwave pulse was 
recorded. It was necessary in such cases to complete 
the groundwave pulse as indicated by the dotted line 
before undertaking the spectral calculations. Since a 
very small fraction of the total area enclosed by the 
waveform is contained under the sketched-in por- 
tion, the error introduced into the spectral calcula- 
tions is small. 

Waveforms of sferics 17 and 18, 3 which were re- 
ceived at ranges of 460 km under nighttime condi- 
tions, are given in figure lb. Note that in these 

3 Numbers were assigned to the atmospherics in order of their selection for 
analyses and each number associates an atmospheric with its spectrum. Sferics 
1, 2, 17, and 18 were selected to represent the variability encountered in the 
reception of the groundwave pulses. The waveforms of the remaining atmos- 
pherics (sferics 3 through 16 and 19 through 33) are not presented in this report. 
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cases the time interval between reception of the 
groundwave and the first-hop skywave pulses is 
about 97 jusec compared to an interval of 147 /xsec in 
figure la. As a consequence of the shorter time 
interval between the reception of the groundwave 
and first-hop skywave pulses, a larger portion of the 
groundwave pulses are concealed by the skywave 
pulses than for the sf erics 1 and 2. Prior to perform- 
ing the spectral analyses of the groundwave pulses 
for sf erics 17 and 18, the trailing edges of the ground- 
wave pulses were sketched in as indicated. A larger 
error resulted in such cases; e.g., the area under the 
sketched-in portion in sferic 18 is approximately 13 
percent of the total area enclosed by the waveform 
of the groundwave pulse; however, tnis error is not 
considered to be serious in the present analyses. 

The ionospheric reflection height for these night- 
time observations was calculated to be about 80 km, 
using a range of 460 km and a time interval between 
reception of groundwave and first-hop skywave 
pulses of 97 /xsec. 

3. Data Analyses 

Let G(t) represent the instantaneous value of tin 4 
vertical component of the electric field in units of 
volts per meter. Also let the waveform start at 
time £=0, so that 

G(t)=0, for £<(), 

and the vwd of the pulse occurs at t=r so that 

G(t)=0, for t>r. 

The complex spectrum of the pulse can be represented 
by the Fourier integral 



FU)< 



£ 



'6(t)dt, 



(1) 



where / represents the frequency in cycles per second 
and t is the time in seconds after the beginning of 
the pulse. 

The complex spectral components of the atmos- 
pheric were calculated at intervals of 1 kc on an 
electronic computer. The absolute value, | F(j) | , of 
the spectra of the groundwave pulses of individual 
atmospherics are given in figures 2a through 2h. 
Fourteen of the analyzed pulses were received under 
daylight conditions from a storm approximately 165 
km from Boulder, and 19 pulses were received under 
nighttime conditions from ranges varying between 
400 and 620 km from Boulder. The approximate 
range at which each pulse was observed is indicated 
on the appropriate figure. 

The frequencies of maximum energy for return 
stroke lightning discharges have been reported by a 
number of different workers.' Their separate findings 
indicate that the frequency peaks occur in the 5- to 
20-kc region. It will be observed that the peak fre- 

4 It is difficult to adequately reference all of the available literature dealing 
with this subject. A few representative references [5-12] have been selected. 
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Figure 2. Amplitude spectra of atmospherics {see text). 



quencies in figure 2 fall within this range of frequen- 
cies. 

With a knowledge of the complex spectrum F(J) 
of a groundwave pulse observed at distance di, and 
of the appropriate values of earth constants, it is 
possible to calculate the spectrum | F(j) | of a ground- 
wave pulse at a different range, d 2) using methods 
evolved by Wait [14]. This rigorous method was not 
followed in this paper, due to the complex nature of 
the calculations and the lack of adequate informa- 
tion concerning earth constants. The less rigorous 
method used in normalizing the spectra \F(J)\ to a 
common distance is described below. 

The atmospheric pulses were recorded at distances 
between 150 and 600 km and their spectral peaks 



201 



were attained at frequencies near 10 kc. A uniform 
earth, having a conductivity of about 10~ 2 mho/m, 
was assumed to exist along the propagation paths. 
Referring to the curves of groundwave amplitude 
versus distance [15], it is evident for these values of 
earth conductivity, distance, and frequency, that the 
signal amplitude decreases approximately at the in- 
verse distance rate. Accordingly, the individual 
spectra were referred to a common distance of 1 km 
by increasing each spectral component in proportion 
to the observation range. The spectrum shown in 
figure 3 is the average of the individual spectra nor- 
malized to this distance. The deviation of the indi- 
vidual spectra from the average is not very large 
which seems to support evidence in the literature 
[12, 13] that the radiation spectra of return-stroke 
lightning discharges do not vary a large amount from 
event to event. The 1-km value was selected as a 
matter of convenience; the spectrum in figure 3 is 
not intended to represent the actual distribution of 
amplitude with frequency at that distance. 
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Figure 3. Average of amplitude spectra, normalized to 1 km. 



The peak amplitude of this spectrum occurs at a 
frequency of 6.5 kc; at frequencies below the peak, 
the spectra] amplitude is proportional to the 0.4 
power of the frequency, while at frequencies higher 
than the peak, the spectral amplitude decreases as 
the inverse 1.5 power of the frequency. 

The total energy available in a lightning discharge 
has been estimated to be about 2X10 17 ergs [16]. 
This large amount of energy is partitioned into a 
number of different forms; for example, some energy 
is required to ionize the gas in the discharge channel, 
a larger portion is liberated as frictional heat in the 
movement of gaseous components, and some is evi- 
dent as electromagnetic radiation, etc. An appreci- 
able part of the energy is consumed in "ground losses" 
due to the flow of induced currents in the earth [17]. 

For the purposes at hand, it was desired to esti- 
mate the amount of vlf radio energy radiated into 
the earth's atmosphere (i.e., exclusive of earth losses) 
by the lightning discharges. This form of energy is 
a very small portion of the total amount available. 



The energy density, E d , of the transient G(t) is 
represented in joules/meter 2 by 



S'-nkf GWdt - 



(2) 



It was assumed that the current at the source was 
a simple linear dipole oriented along the vertical axis. 
The radiation component of the dipole electric field 
is proportional to cos 0, where 6 is the angle of eleva- 
tion measured from the horizontal. The vlf energy, 
E T , in joules, radiated into the atmosphere was ob- 
tained by integrating the energy density E d over the 
surface of the radiation pattern of the source. Thus, 



E. 



Jo Jl 



7T/2 



E d r 2 cos 3 dddd<t> 



(3) 



=4l3wr 2 E d , 



where 8 is the angle of elevation measured from the 
horizontal, <j> is the angle measured in the equatorial 
plane, and r is the distance in meters from the source 
to the point of observation. 

Values of energy represented by each of the pulses, 
calculated using expression (3), are given in table 1. 

Table 1. Very-low-frequency energy radiated into the atmos- 
phere by return-stroke lightning discharges 



Atmos- 
pheric No. 




E T 



Sferic groundwave recorded at Boulder 
during daylight July 7, 1959 



1. 
2_ 
3. 
4_ 
5. 

6. 

7. 
8. 
9. 

in 

11 
VI 
L3 

11 



km 


kc 


165 


10.0 


165 


5.0 


165 


9.0 


165 


7.0 


165 


9.0 


165 


5.5 


165 


9.5 


165 


10.0 


165 


6.0 


165 


6.0 


165 


7.5 


165 


9.5 


165 


10.0 


165 


9.5 



2.2X19 3 
6.2X103 
1.7X10* 
5.5X10* 
1.8X103 

2. 1X10* 
1.1X102 
1.5X10* 
4.3X103 

8.8X10* 

6.9X102 
1.3X105 
1.5X103 
5.9X102 



Sferic groundwave recorded at Boulder 
during nighttime July 26, 1958 



L5 
16 
31 
33 

17. 
18 
19 
20 
2] 
32 

22 
23 

24 
25 
26 

27 
28 
29 

30 



400 


7.5 


400 


6.0 


400 


8.5 


400 


8.5 


460 


6.5 


460 


5.5 


460 


5.5 


460 


7.0 


460 


7.0 


460 


5.5 


600 


15.5 


600 


12.0 


620 


5.5 


620 


7.0 


620 


6.5 


620 


9.0 


620 


8.5 


620 


8.0 


620 


6.5 



2.1X10* 
6.5X103 
2.0X10* 
5.5X103 

4.5X103 
2.9X10* 
6.2X103 
3.1X10* 
1.1X10* 
9.3X103 

7.9X103 
1.1X10* 

1.1X105 
9.2X10* 
6.4X10* 
1.3X10* 
9.6X103 
6.4X103 
6.3X10* 
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The average of these values of energy is 26,600 j. 
On another occasion [10], sferics were received from 
thunderstorms in a frontal region extending through 
Texas, Louisiana, and Oklahoma. The average of 
the total radiated energy values calculated from 
seven of these sferics was about 250,000 j ; the indi- 
vidual values range from 130,000 to 500,000 j. Laby, 
et al. [11], reported the average value of energy radi- 
ated from 35 sferics observed in Australia was 200,000 
j. It appears that the sferics used in the present 
work resulted from less energetic discharges than 
those reported by Taylor and by Laby. It is not 
at all certain that the energy value reported here is 
representative of return-stroke discharges from thun- 
derstorms occurring in the Rocky Mountain area. 
However, part of the differences between the energy 
values quoted might have resulted from differences 
in thunderstorm characteristics, differences in topog- 
raphy, earth conductivity, etc. 

The importance of land elevation in determining 
thunderstorm characteristics was pointed out by 
Loeb [18]. In general, the lightning discharges used 
in this analysis occurred over land of higher eleva- 
tion than those previously reported by Taylor. It 
is not known in what area of Australia the lightning 
flashes occurred that were used in Lain 's analysis. 

It is interesting to note that the ratio of the total 
energy available in a lighting discharge to the values 
of electromagnetic energy radiated into the atmos- 
phere is as large as 10 5 or 10 fi . 

A wide variation was noted in the detail of the 
individual ground wave pulses; however, most of them 
consisted of two half-cycles. The first half-cycles 
were larger in amplitude and shorter in duration 
than the corresponding values in the second half- 
cycle. These two characteristics of the first half- 
cycle would be expected to control the total radiated 
energy to a certain extent. After normalizing the 
observations to a 1-kni range, the relationships 
existing between various parameters mentioned above 
were examined. Some of these comparisons are 
described below. 

(A) Frequency of spectral peak versus total radiated 
energy. The relationship between the frequency of 
spectral peak (/ ) and the total radiated energy (E T ) 
is given in figure 4. Although the data points 
scatter, there is a slight tendency for the more 
energetic sferics to attain energy peaks at the lower 
frequencies. 

(B) Peak amplitude oj sj erics versus total radiated 
energy. The peak amplitude of the groundwave 
pulse (E m&x ) is plotted as a function of total radiated 
energy in figure 5. There is a close relationship be- 
tween these parameters. The straight line drawn 
through the data points relate total radiated energy 
to the 2.4 power of the peak field strength. 

(C) Spectral peak versus total radiated energy. The 
relationship between the peak of the spectrum of 
the groundwave pulse, E(f) m&x , versus the total radi- 
ated energy of the pulse is given in figure 6. A close 
relationship is indicated between these parameters. 
The straight line through the data points indicates 
that the energy is proportional to the 1.8 power of 
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Figure 4. Frequency of spectral peak versus total radiated 
energy. 



I0 J 


















«*"^ 










J> 


o ~,^^° 






10* 






^^^ 










^^ 















E T = 0.012 E 2 m 4 Q 2 x 








in 














,o 2 


io 3 


o 4 IO 5 


,o 6 



E T (JOULES) 

Figure 5. Pulse amplitude versus total radiated energy. 
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Figure 6. Spectral peak versus total radiated energy. 

the spectral peak. The frequency of the spectral 
peak is not explicitly accounted for in (his relation- 
ship. 

(D) Length of first half -cycle versus total radiated 
energy. In figure 7 the relationship between the 
length of the first half-cycle of the groundwave pulse 
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Figure 7. Length of first half-cycle versus total radiated energy. 

(t x in microseconds) and the corresponding value 
of total radiated energy is plotted. The relative 
importance of the first half-cycle in determining 
the energy is indicated in this plot and in figure 5, 
in which the peak amplitude of the groundwave 
pulse is plotted as a function of total radiated 
energy. 

The field strength of a pulse is related to the rate 
of change of the source moment [13, 14, 19], hence, 
flashes having rapidly rising currents would be 
expected to produce sf erics of high peak amplitude. 
The duration of the first half-cycle of the ground- 
wave pulse should correspond to the time at which 
the source current reached its maximum value. 
Hence, sf erics of large energy result from rapidly 
rising and long continuing source currents. It 
would be expected that the more energetic discharges 
would produce sf erics of longer duration, and a 
strong tendency of this nature is observed in the 
data; however, it was also observed that sf erics 
representing relatively low values of energy attain 
spectral peaks within the same frequency region, 
5-20 kc, as those of larger energy. 

The close relationship existing in figure 5 between 
the peak amplitude of the first half-cycle and total 
radiated energy suggests the possibility of gaining a 
useful estimate of energy by observing the peak 
amplitude of the groundwave pulse and by de- 
termining the distance to the flash. However, it is 
not known that this relationship will hold equally 
well at the higher energy levels. 

It would be of interest to compare the character- 
istics of sf erics originating (1) over land of different 
elevations and those originating over sea water, 
and (2) from frontal and air-mass thunderstorms. 
This would aid in appraising the influence of topog- 
raphy, earth conductivity, and meteorological factors 
in determining the characteristics of sf erics. 

Atmospheric waveforms recorded at Boulder were 
used in the present analysis. Simultaneous record- 
ings of these atmospherics are also available from the 
other stations located at Stanford and Salt Lake 
City. Analyses of this data is continuing to de- 



termine values of propagation attenuation [10] and 
phase velocity of propagation [20, 21]. 
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